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Abstract 
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theories based on o[^^ have real spectra in finite volume, the spectra of a^2^ models is 
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1 Introduction 



Unitarity plays an important role in quantum field theory (QFT). In many application of 
QFT, the framework of quantum theory requires a positive definite conserved probability, 
which is guaranteed by the Hermiticity of the Hamiltonian and the unitarity of the S 
matrix. Hermiticity also guarantees that the spectrum of the Hamiltonian is real, an 
important condition if the Hamiltonian is to be interpreted as the energy of the physical 
system. 

However, there are applications of the QFT formalism when this is not a necessary physical 
requirement: non-unitary QFTs (even those with complex spectra) appear to play an in- 
creasing role in the investigation of statistical mechanical systems (e.g. disordered systems 

i)- 

In this paper we intend to study the simplest cases: the finite-size spectra of models 
related to imaginary coupled Toda theories a^^^ (sine-Gordon) and a^2^ (ZMS), with periodic 
boundary conditions. Our main aim is to establish conditions under which these theories 
have real spectra, continuing our earlier work started in |^. Reality of the spectrum 
could be interesting for several reasons: (i) it affects the large distance asymptotics of 
correlation functions (in case of a complex spectrum, these asymptotics can show oscillating 
behaviour), (ii) it could allow a redefinition of the theory that renders it unitary (although, 
as we discuss, physical requirements may prevent that, and it is not at all clear whether this 
can be performed while maintaining the interpretation of the given model as a local QFT). 
Finite size spectra are also interesting in their own right: periodic boundary conditions 
can also be thought of as realization of finite temperature. 

We start by recalling the general issues involved, and the relationships between the three 
distinct concepts of probability conservation, unitarity and the reality of the spectrum. 

For any imaginary coupled affine Toda field theory (ATFT), there are three different classes 
of models which can be considered, namely the original (unfolded, unrestricted) , the folded 
and the restricted (RSOS) models. We define these in section ^ and discuss the reality of 
the spectrum in these three classes for the theories related to a^^^ (sine-Gordon) and a^2^ 
(ZMS) in turn. 

In the case of a^^\ the finite-size spectrum of the RSOS models turns out to be a subset 
of the spectrum of a suitable folded model - i.e. the state space of the RSOS model can 
be thought of as a projection of the space of the folded model. Using this relation, we can 
conclude that the spectra of all a^i^ related models are real. In this paper we also give 
evidence that a similar relationship exists between the RSOS and folded models based on 
02 . However, this is not sufficient to establish the reality of these models. We then turn 
to transfer matrix arguments and numerical analysis to investigate the spectrum. 

In all cases in which our studies find no violation of reality, there are several possibilities 
remaining: 



2 



1. Apart from cases in which the S'-matrix can be shown to be Hermitian analytic (e.g. 
perturbations of unitary minimal models), we only have numerical (and sometimes 
matrix perturbation theory) results. Since these numerical tests require 'scanning' 
over the full range of rapidity of each particle in a multi-particle state, there is always 
a chance that such methods miss some region of rapidities in which reality is violated. 
In addition, numerical diagonalisation always introduces some numerical errors, and 
so some threshold condition must be defined to separate real cases from nonreal cases. 
It remains a question whether reality violations below the threshold are genuine or 
not. 

2. We could only examine transfer matrices containing up to five particles, as the nu- 
merical calculations become progressively more and more difficult as we increase the 
number of particles involved. Therefore, it is possible that higher particle transfer 
matrices would introduce further constraints on the reality of the spectrum. 

(2) 

3. In the case of 02 , our methods only allow us to determine the large volume asymp- 
totics of the spectrum which have a power-like decay as a function of the volume. It 
is possible that further contributions (decaying exponentially with the volume, e.g. 
those related to vacuum polarisation) would spoil reality. 

(2) 

4. For 02 , the bootstrap is closed only in certain regimes of the parameter space and 
we do not know the S'-matrices of all possible particles in the spectrum. It is possible 
that reality of the spectrum is only violated in sectors which contain such 'extra' 
particles. 

As a result, we can only have definite results in the cases when the spectrum is complex: 
if we find that some multi-particle transfer matrix has non-phase eigenvalues then we can 
conclude that the Hamiltonian must have complex eigenvalues. Whenever our examination 
finds that the spectrum is real we cannot the possibility exclude that further study would 
find a complex spectrum. 

However, the theories for which we cannot find any violation of reality show some very 
striking patterns, and therefore we believe it is possible that several or possibly all theories 
falling into these patterns have real spectra. Clearly, further understanding is necessary: 
the most promising approach would be the development of some exact method to describe 
the finite size spectra (e.g. an extension of the 03 NLIE [|] to describe excited states) and 
a deeper understanding of the relation between the folded and the RSOS models in finite 
volume. 

As a side result, we also correct some minor mistakes in our previous paper and the 0i5 
S'-matrices conjectured earlier by one of us in 0|. We summarise our results in section ^. 
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2 Unitarity and related concepts in QFT 



In this section we discuss the relation between the concepts of probability conservation, 
unitarity and reality of the spectrum in quantum field theory which are often confused. 
We also recall the definitions of Hermitian analyticity and R-matrix unitarity and their 
relationship with unitarity of the S'-matrix. 

2.1 Unitarity, reality and probabilistic interpretation 

For all the theories we consider there exists a non-degenerate sesquilinear form (which we 
call an inner product for short, even if it is not positive definite) on their space of states 
which is conserved under the time evolution described by the Hamiltonian i.e. the time 
evolution operator and the S'-matrix preserves this inner product (and the Hamiltonian is 
Hermitian with respect to the conjugation defined by it). For perturbed conformal field 
theories such an inner product is inherited from the standard inner product used in CFT, 
even away from the critical point. 

In the cases when this inner product is positive definite (and defines a Hilbert space struc- 
ture on the space of states) this implies the usual Hermiticity/unitarity, with the conse- 
quence that the S'-matrix has phase eigenvalues and the energy spectrum is real. This inner 
product makes possible the usual probabilistic interpretation in quantum theory. These 
theories are called unitary QFTs. 

On the other hand, in many systems the inner product is indefinite. It is still possible for 
the spectrum to be real and the S-matrix eigenvalues to be phases in which case we call 
the theory a real QFTs. However, generally when the inner product is indefinite, then the 
spectrum is complex and the eigenvalues of the S'-matrix are not phases, in which case we 
call the QFT nonrea^. 

One can see that the three notions of unitarity, reality and probability conservation are 
in fact different. If we want the theory to have a positive definite conserved probability, 
then unitarity is necessary, which is a stronger notion than reality. However, often this 
is not required; moreover, there exist exotic probability theories which allow for negative 
or complex probabilities (which in that case of course loses the interpretation in terms of 
frequency of events). 

An important point is whether the S'-matrix has any relevance in the case when the theory 
is non-unitary. It is easier to accept this when the theory is still real, as in the case of the 
scaling Lee- Yang model (i.e. Virasoro minimal model (2, 5) perturbed by its single non- 
trivial relevant operator (pi^) or for indeed of any 0i3 perturbations of Virasoro minimal 
models. However, there are examples of nonreal theories (e.g. the Virasoro minimal model 

^ There doesn't seem to be any agreed convention for the naming of what we call real and nonreal 
theories. We simply use these names in this paper for convenience. 
^For an early discussion of this issue, cf. 



4 



(3, 14) perturbed by operator 0i5 0) for which the finite size spectrum extracted using 
thermodynamic Bethe Ansatz (TBA) for the vacuum energy and the Bethe-Yang equa- 
tions for multi-particle states matches perfectly the spectrum of the Hamiltonian obtained 
numerically using truncated conformal space approach (TCSA), even for the complex part 
of the spectrum. 

To summarise, we believe that for a large class of theories (unitary, real and non-real), and 
in particular for the RSOS models considered in this paper, the spectrum is determined 
by the S-matrix even when the spectrum is complex. 

2.2 Unitarity and Hermitian analyticity 

In analytic S-matrix theory the property of unitarity is closely linked with that of Hermitian 
analyticity ||^, ||]. Without entering into details, we intend to recall these concepts and their 
relations here as they are going to play an important role later. 

Unitarity is simply the statement that using a Hermitian conjugation f with respect to a 
positive definite inner product, the S operator that maps out-states into the in-states has 
the property 

SS^ = S^S = I. (2.1) 

Given a theory that is non-unitary but real, since the S operator has phase eigenvalues 
it is obviously possible to define a new positive definite inner product with respect to 
which S is unitary. However, this inner product may be inconsistent with the rules of 
analytic S matrix theory and/or may render the Hamiltonian non-Hermitian (while the 
Hamiltonian was Hermitian with respect to the original, indefinite inner product). The 
scaling Lee- Yang model is an example of this situation. Its spectrum contains a single scalar 
particle which appears as a bound state in the two-particle scattering. S'-matrix theory 
then relates the norm of one-particle states to two-particle states through the residue 
of the corresponding pole in the two-particle S-matrix. The natural inner product in 
the perturbed CFT formalism (defined through the Hermiticity of Virasoro generators) 
is indefinite, with n-particle states having the signature (—1)" (or —(—1)", depending on 
choice, in which case the sign of the residue is consistent with S'-matrix theory. However, if 
one attempts to use a positive definite inner product, then the residue of the two-particle 
S-matrix at the bound state pole has the "wrong sign", which is how it is often stated in 
the literature. As a result, the natural inner product of the scaling Lee- Yang model is 
indefinite, and it is an example of a real but non-unitary theory. Therefore the existence 
of a positive definite redefinition of the inner product does not mean the theory can be 
made unitary because it may confiict with some other physically motivated requirements. 
For a further discussion of this issue see [|]. The property that the residues of the bound 
state poles have the "right sign" is sometimes called "one-particle unitarity" 0. 

In integrable theories the whole S operator is encoded in the set of two-particle S-matrices 
Sab{(^), where A and B denote the particles (or multiplets if there are internal quantum 
numbers) and 6 is their relative rapidity. It is a simple matter to prove that unitarity of 
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all possible two-particle 5'-matrices 



'^SABid)ij (5'ab(6')^„) = SrrnS^ 



(2.2) 



(where we explicitly wrote the multiplet indices) implies the unitarity of multi-particle 
transfer matrices (defined in Appendix 0), and also means that the spectrum is real. In 
writing equation (|2.2| ) we assumed that we had chosen an orthonormal basis in the internal 
multiplet space. 

If, in addition, the inner product on the space of states is positive definite then the theory 
in question is a unitary QFT. We shall abbreviate the property ( p. 21 ) as TU (two-particle 
unitarity) . 

Hermitian analyticity (HA) tells us something about the behaviour of the S'-matrix ele- 
ments under complex conjugation. It states that 



On the other hand, the S'-matrices we investigate here are derived from quantum group R- 
matrices. -R-matrices also satisfy a relation known as "unitarity" in quantum group theory 
which we call here i?-matrix unitarity (RU) and takes the form 



We see that RU and HA together imply TU and thus reality of the spectrum (note that 
equation (|2.2|) is meant only for physical i.e. real values of 9). Therefore for the models we 
consider Hermitian analyticity implies a real spectrum since RU is automatically satisfied 
due to the quantum group symmetry. 

These notions can be appropriately generalised to RSOS theories, where the multi-particle 
polarisation spaces are not simply tensor products of one-particle ones. Rather, the multi- 
plet structure is labelled by so-called RSOS sequences, which denote 'the vacua' between 
which the particles mediate and are constrained according to so-called 'adjacency rules'. 
The two-particle S'-matrix therefore carries four vacuum indices a, b, c, d and it is specified 
in the following way (we omit particle species labels for simplicity): 



(2.3) 




(2.4) 



e 




(2.5) 
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For RSOS theories TU, HA and RU take the form: 



TU: 5^5^^(0)5^f(0)* = 5,, (2.6) 

d 

HA : SfM* = STA-e*) (2.7) 
RU: Y.^^MSTd{-0) = Sbe (2.8) 

d 

3 The models 

The models we study can be classified in the following way: 

1. The original (unfolded, unrestricted) models. These have an infinite set of degenerate 
'vacua' and their spectra are built from a fundamental soliton doublet (a^i^) or triplet 
{a'^^) by closing the S-matrix bootstrap. 

2. Folded models. Using the periodicity of the field theoretic potential, one can choose 



to identify the ground states after k periods (see ||T0| for the sine-Gordon case), i.e. in 
a fc-folded model one has k ground states, between which the solitons mediate. The 
spectrum and the scattering theory for this case can be straightforwardly written 
down using the well-known S'-matrices of the original model. 

Restricted (RSOS) models. At certain 'rational' values of the coupling it is possible 
to define a space of 'RSOS states' as a quotient of the full space of states on which 
the action of the S operator is well defined. These states can be labelled by sequences 
of 'vacua' and the S'-matrix factories into 2-particle 'RSOS type' S'-matrices. These 
RSOS S'-matrices describe the scattering theory of the 0i3 or 0i2/02i/0i5 perturba- 
tions of Virasoro minimal models, respectively for a^i^ and a^2^. These restrictions 



were discussed in the following papers: the 0i3 case in ||TI|, 0i2/02i in ||T2| and 0i5 in 
[|]. The above construction of the RSOS states relies on the action of the quantum 
group and has only been defined for the theory on the full line i.e. in infinite spatial 
volume. However, given the two-particle S'-matrix derived in this way one can easily 
define the theory in finite volume. 

All of the models enumerated above are integrable and therefore there are numerous results 
for their finite volume spectra. We shall mainly use results obtained using two methods. 

Firstly, the so-called nonlinear integral equation (NLIE) approach, which gives exact results 
for the spectrum and was pioneered by Kliimper et al. ||T3| on the one hand and Destri 
and de Vega fl^ on the other. The original developments concerned mainly the spectra of 
the a^^'* case (for the excited states see also ||15], [T^). 
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Secondly, we use the related method of the Bethe-Yang equations (cf. and references 
therein), which describe the large volume asymptotics of the finite size spectrum. If the 
NLIE for a given theory is known, then the asymptotics can be derived independently and 
for a^i^ it was checked that they agree with the results from the Bethe-Yang equations 



T7[| . In general when the NLIE is not known this is the only known analytic way 
to obtain information about the finite size spectrum (in some cases, TBA equations are 
known for the excited states [|l^], but these cases are covered by the NLIE as well). 

In finite volume, since the field theory potential is periodic, it is possible to introduce 
twisted sectors of the original model ||20|- These are labelled by a twist parameter a, 
defined mod27r. The twisted sectors have different finite size spectra which are given by 



a modification of the original NLIE in which a appears explicitly |^ . 

The finite size spectrum of the fc-folded model is the union of the spectra of twisted sectors 
where the twist runs through the values |[T0| 



27rm 

a = — — , m = 0,...,k — l (3.1) 

fV 

When we turn to specific models we shall discuss a relation between the finite size spectra 
of folded and RSOS models, namely that the finite-size spectrum of an RSOS model is 
exactly a subset of the spectrum of an appropriate folded model. Such a relation first 
emerged for the case of a^^^ 0, and we shall give evidence that a similar relationship 
holds for a'^\ 

To be precise, for a^^-* it is known that both the exact spectra (described by the NLIE) 
and, as a consequence, their asymptotics (described by the Bethe-Yang eqns) of folded and 
RSOS models are related. For 03 it is known that the exact vacuum energies in finite 
volume are related (this was shown using the NLIE framework in 0]), but at present the 
NLIE for excited states is not known. We shall present evidence that the asymptotics of 
the excited state energies given by the Bethe-Yang eqns are similarly related. We believe 
similar relations will hold for the exact spectra of every ATFT and the corresponding RSOS 
models. 

There exists some evidence in support of this claim. Al.B. Zamolodchikov has shown that 
within the perturbed conformal field theory framework, the perturbative expansion for 
the vacuum energy of the sine-Gordon model with a suitable value of the twist parameter 
could be reinterpreted as a perturbative expansion in the RSOS model |^], to all orders 



in perturbation theory. There is also a quantum group argument for the agreement of the 
vacuum energies ||22|, which relies on a formulation of the partition function in terms of 



the system in infinite volume. In addition, using the a^^^ NLIE it is possible to calculate 
(a) the exact ultraviolet conformal weights and (b) the energies of the ground state and 
excited states in the twisted sectors of sine-Gordon model to very high numerical accuracy. 
Comparison of these data to (a) the known CFT data and (b) numerical finite volume 



spectra extracted using TCSA, respectively, shows an excellent agreement |21 . 
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It should be possible to find a proof using quantum group arguments directly for the theory 
in finite volume, by finding an appropriate projection from the folded model to the RSOS 
space which commutes with the transfer matrices, but we are unaware of such an argument. 
The knowledge of this projection would also be useful because it could be used to select 
systematically the RSOS states in the NLIE approach. 

We now proceed to give our conventions for the models based on a^^^ and a^'^ and give the 
RSOS — folded relations that we shall check in Section ^. 



3.1 Conventions for a\ 



(1) 



The a^i'' theory is simply the sine-Gordon model, for which we take the action to be 



A.. 



sG 



X 



1 ?77^ \ 

-a,$a^$ + -^cos/5<l>j . 



(3.2) 



As is well known, introducing the parameter q = exp^Sni / (3^) ^ one can show that the model 
is invariant with respect to Uq{a^i^). The spectrum consists of a doublet of solitons and a 
collection of scalar particles (breathers). As described by the RSOS theory A^py + 0i3 
is obtained as a restriction at /? 

m 



the results of 17 and 



Sirp/p', that is at g = exp{7rip' / p) . Putting together 
it is straightforward to see that the finite volume spectrum of 



this RSOS model is a subset of that of the 2p-folded model. 



3.2 Conventions for a 



(2) 



(2) 

We take the action of the a\ model to be 



A 



ZMS 



d'x 



1 TTl^ / \ 

+ ^ exp (2iy^$) + 2 exp {-iy/j<!>) J 



Introducing the parameter 



q = exp 



(3.3) 
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(3.4) 



one can show that the model has a symmetry under the affine quantum group Uq{a2 ) and 



as a result its S'-matrix can be explicitly constructed [12 



The spectrum of the original model consists of a fundamental soliton transforming as a 
triplet under Uq{a^i^), having an S'-matrix of the form 



Sqq{x,'j) R{x, q) , where x = exp{2n9 / ^) and ^ 



2-E 7 



(3.5) 



3 (27r-7) ' 

where S'qq is a scalar function which is always a pure phase for real x and 7, and R{x, q) is 

(2) 

the i?-matrix of Uq{a2 ) in the triplet representation. 
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Depending on 7, there may be other particles in the spectrum. For ^ < vr there are 
breathers Bq , n = 1, . . . , [tt/^] formed as bound states of Kq transforming in the singlet 
representation, and for ^ < 27r/3 there are higher kinks Ki , n = 1, . . . , [27r/(3^)], also 
formed as bound states of Kq and transforming in the triplet representation, and their 
associated breathers B^. 

In addition, it is known that there are values of 7 for which the bootstrap does not close 
on this particle content, and for which there are more particles in the spectrum. Since we 
are only able to investigate the transfer matrices for the kinks Ki and breathers S", it is 
only possible to prove that any particular model has a complex spectrum. If we find that 
the spectrum for these particles is real, that does not imply the reality of the full spectrum 
simply because the S'-matrix bootstrap could still generate S'-matrices with non-phase 
eigenvalues. 

The S'-matrices of the higher kinks have the form 

= S^,„(x,7)i?((-ir+"x,g) (3.6) 

where is a scalar factor. 

Since the scalar factors Sjjj „ are phases and there are only two different matrix structures 
for the two-particle S-matrices SK,ri,K„ i-e R{x ,q) if m + is even and R{—x ,q) if m + n 
is odd, for the purposes of finding out if the eigenvalues of the transfer matrices (defined 
in appendix [C]) it is hence to sufficient to consider 

• For Tc/j < 1: transfer matrices containing only Kq 

• For 7r/7 > 1: transfer matrices containing only Kq and/or Ki 



3.2.1 RSOS models related to a 



(2) 



The a^2^ model has two inequivalent RSOS restrictions ||T2|, ||. From the conjectured ground 
state NLIE of the a2^''-related models § one can determine the minimal value of the folding 
number for which the RSOS ground state is in the spectrum of the folded model. Together 
with the results in 



12, HI, it leads to the following conjectures 



1. Aip^p' + 012 is a projection of the p-folded ZMS model at 7 = irp/p' . 

The RSOS S'-matrices are given in the appendix |^ where we have specified all the 
formulae including the 6j symbols since the ones we found in the literature all had 
typos.0 

As for the unrestricted case, the matrix part of the S-matrices depend only on q so 
they depend on p' only modulo 2p. In fact, we only need study the values 1 < p' < p 
since the substitution p' — > — p' changes the S-matrix elements to their complex 

^We are grateful to Giuseppe Mussardo (SISSA) for providing us with a correct set of formulae. 
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conjugate (if we choose the square root branches in the 6j symbols appropriately) 
and thus simply complex conjugates the transfer matrices and their eigenvalues too 
as well. 

(n.b. the models A^py + ^21 are contained in this class through the identification 

Mp^p' + (f)2l = Mp'^p + 012.) 

2. A4p^p' + 015 is a projection of the 2p-folded ZMS model at •y = Anp/p'. 

The RSOS restriction leading to the S'-matrices was performed in f^, however, certain 
amplitudes had wrong prefactors. The reason was that in the scattering amplitudes 
of two (oppositely) charged solitons into two neutral ones there is a Clebsh-Gordan 
factor which was not taken properly into account and therefore the amplitudes of 
[|] do not satisfy the Yang-Baxter equation and i?-matrix unitarity. The corrected 
amplitudes are listed in Appendix |B[ 

These perturbations lead to renormalisable field theories if the conformal weight of the 
perturbing field is less than or equal to two, which is equivalent to the condition n/j > 1/2 
for both cases, which is also the condition for the unrestricted model to be well-defined. 

4 Results 

The leading approximation (e.g. in the NLIE formalism) to finite size effects in large volume 
is given by the Bethe-Yang equations which are summarised in Appendix y. In particular 
it can be seen that the spectra of the models can only be real if the transfer matrix 
eigenvalues (denoted in the Appendix by A*^*-* . . . ,'0^)) are all phases for real values 

of the particle rapidities 

We check for which values of the coupling the various transfer matrices have phase eigen- 
values. There are some simplifications we can make. 

1. In all our calculation we omit the scalar prefactors of the S'-matrices since these are 
irrelevant for determining whether the eigenvalues are pure phases or not. We call the 
transfer matrices obtained from the S'-matrices without the scalar factor 'reduced' 
transfer matrices. 

2. The i?-matrices of all models we consider have some discrete symmetries (which we 
describe here for real values of the rapidity) . 

(a) R{x, q*) = R{x, g)*, which means that we only need to consider < arg(g) < vr. 

(b) R{x, —l/q) = UR{x, q)U^^ where t/ is a diagonal matrix whose nonzero entries 
are ±1. 
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The reduced transfer matrices are still in general very complicated, even for the original 
models, and therefore in our study we diagonalised them numerically for a large set of values 
for the rapidities and other parameters. In all our plots we show only the eigenvalues of 
the reduced transfer matrices. 



4.1 a^^^ related models 

In the case of a^^^ related models, since the folded and original models are all unitary, their 
finite size spectra are real and therefore all 0i3 perturbations have real spectra as well, 
regardless of whether they are perturbations of unitary or non-unitary minimal models. In 
addition, the NLIE formalism ||21| yields manifestly real spectra as well. 

In other words: 

., r,,, ., RSOS real 

origmal unitary ^ folded unitary < , ^r^^^^ 

^ i 7^ RSOS unitary 

One might think that, since the RSOS spectrum is simply a subset of the spectrum of 
the unitary unrestricted model, the RSOS model would simply inherit a positive definite 
inner product and also be unitary. However, the fact that the spectrum of the RSOS and 
unrestricted models are different, means that the constraints of analytic S-matrix theory 
may enforce different inner products. As an example, consider the value Sn/P^ = 2/5. The 
particle content of the unrestricted model is a soliton doublet s, s and a single breather 
B. The RSOS model, on the contrary, is the Lee- Yang model with a single particle B 
(the solitons are removed from the spectrum by the RSOS restriction). There is a first 
order pole in Sbb{0) at ^ = 27ri/3 which must have an explanation in terms of on-shell 
diagrams. In the unrestricted sine-Gordon model, it is explained by the famous Coleman- 



Thun mechanism ||23|| as the sum of singular contributions from diagrams with internal 
soliton lines. In the RSOS model there are no solitons so this pole must be explained by 
a single diagram in which B occurs as a bound state of two B particles. As we discussed 
before, the sign of the residue of this pole forces the inner product to be indefinite. 

By numerical diagonalisation of the reduced transfer matrices we found that the transfer 
matrix eigenvalues for the RSOS models A^^y + 0i3 are a subset of those for the 2p-folded 
sine-Gordon model with /3 = \f%Tvp~fp. We also observed that the numerically computed 
transfer matrix eigenvalues are all phases, as expected from the general argument above. 

In figure |l], we plot the arguments of the eigenvalues of the reduced two-particle transfer 
matrices of the RSOS models (A^7,m + ^is) and the 14-folded sine-Gordon model. The 
RSOS models are only defined for m integer, i.e. for arg(g)/7r taking values m/7. For 
these values the eigenvalues of the transfer matrices (shown as blobs) are a subset of the 
eigenvalues of the matrices of the folded model shown as lines. 

The two-particle space of the 14-folded model has dimension 28, and the spectrum consists 
of 12 eigenvalues each of multiplicity 2 and 4 of multiplicity 1. The RSOS model excludes 
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Figure 1: Comparison of the eigenvalues of the matrix part of the two-particle transfer matrices 
for the RSOS model {Mj^m + '/'is) and the 14-folded af ^ ATFT The horizontal axis is arg(g)/7r 
restricted to the fundamental domain between and 1, and the vertical axis is the argument of 
the eigenvalues. The folded models are shown in green/solid lines, and the RSOS models (which 
are only defined for discrete values of q) are shown as blobs. The relative rapidity of the particles 
was chosen = 1/9. 



the eigenvalues of multiplicity 1, 2 of the other eigenvalues, and each of the remaining 10 
eigenvalues appears with multiplicity 1 only. 

There are RSOS models A^py + 0i3 for which the 5'-matrix is Hermitian analytic (e.g. the 
unitary case p' = p + 1 or the series p = 2, see ||ll|) which then implies that the spectrum 
is real. Our arguments above, however, extend to all 0i3 perturbations. 

(2) 

4.2 a2 related models 

(2) 

In this section, we investigate models related to . We start with the 1-folded model 
where we correct a mistake in We then consider multi-particle transfer matrices and 
show the emergence of a pattern which seems to carry over to the folded case and it plays 
an important role in the RSOS well. 

For all three classes of model, we must consider separately the case 71/7 > 1 for which 
there are higher kinks, and 71/7 < 1 for which there are no higher kinks and consequently 
fewer constraints on the values of 7 allowed by the reality of the spectrum. 

We recall that due to the symmetries of the i?-matrix, the spectrum of models with 
— TT < arg(g) < is obtained by conjugating the spectrum of the model with 1/g, and so 
we can restrict our attention to models with < arg(g) < tt. 
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(2) 

4.2.1 The 1-folded models based on 

We shall consider first the two-particle transfer matrix, and then the higher particle number 
transfer matrices in turn. As we have mentioned, for our investigation of the reality of the 
spectrum, we can ignore the scalar prefactors in the S'-matrices and need only consider 
the 'reduced' transfer matrices constructed out of the appropriate /^-matrices. The two- 
particle 'reduced transfer matrices' are simply the i?-matrices themselves: if the kinks 
being scattered are Km and i^n, then the appropriate i?-matrix is l)™+"x, g). 

As calculated in 0, the eigenvalues of R{x,q) are three pairs of doubly degenerate eigen- 
values, 

'l-q^Vx\ [l + q^^\ ^^^^ 



and three eigenvalues A satisfying 



(x-g^)(a; + g6)A3 + q^{2 \ q^){\^ \ x'^X) 

+ x(g2-l)(l-3g4 + g8)(A2 + A) (4.2) 
- g2(i + 2g2)(x2A2 + A) + (1 - g^x)(l + g^x) =0. 



As stated in 0, for a; negative, there are obviously eigenvalues which are not phases, 
except when = 1. However, in [|] we stated incorrectly that all the eigenvalues were 
phases for x non-negative and q a phase, and instead the correct result is that there are 
non-phase eigenvalues for < | arg(g)/7r| < 1/4. This is due to branch cuts that appear in 
the solutions of the cubic equation ( |4.2| ) which we overlooked before. 

We were not able to diagonalise any higher particle number transfer matrices exactly. 
Therefore we resorted to numerical methods for the cases of three, four, five and six par- 
ticles, combined with matrix perturbation theory in the case of three and four particles. 

The 1-folded models with 7r/7 > 1 

If 7r/7 > 1 then the model contains both and K^. Since the matrix part of Sko,Ki is 
R{—x, g), the eigenvalues are just given by (|4.1|) and (|4.2|) with x replaced by — x, and from 
looking at (|4.1| ) we see that the two-particle transfer matrix has non-phase eigenvalues 
unless g^ = 1. This means that, unless g'^ = 1, the theory is nonreal. We have not 
investigated the case g^ = 1 any further. 

The 1-folded models with n/j < 1 

If 7r/7 < 1 then the only particles in the spectrum are the fundamental kinks (and the 
first breather for 7r/7 > 2/3 ). In particular, the model contains no higher kinks, and the 
only reduced transfer matrices to be considered are those constructed out of R{x, q) with 
X positive. 
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The eigenvalues of the two-particle case have already been presented, and show that the 
theory is nonreal for < |arg(g)/7r| < 1/4 and for 3/4 < |arg(g)/7r| < 1. 

Numerical investigation of the three-particle transfer matrices show that the theory is 
also nonreal if 1/4 < |arg(g)/7r| < 1/3 or 2/3 < |arg(g)/7r| < 3/4. It turns out that 
the strongest violation of reality happens for small x and and one can expand the 
eigenvalues in a perturbation series around x = ?/ = 0. Some care must be taken as 
the X = y = Q transfer matrices have a nontrivial Jordan form, so the Lidskii-Vishik- 
Lyusternik generalised perturbation theory must be used for the expansion. Without 
entering into technicalities we only wish to mention that one obtains exactly the same 
pattern as described above. 

Consideration of the four-particle transfer matrix adds further regions of non-reality, 1/3 < 
|arg(g)/7r| < 3/8 and 5/8 < |arg(g)/7r| < 2/3 , and of the five-particle transfer matrix 
additional regions of 3/8 < |arg(g)/7r| < 2/5 and 3/5 < |arg(g)/7r| < 5/8. 

Assuming this pattern to continue for higher particle number transfer matrices, we are led 
to the conjecture that the eigenvalues of the n particle transfer matrices are always phases 
for 

arg(g) 



1 

< 

2n - 



vr 



and at the isolated points 



arg(g) 



TT 



2 2m 



1 1 
< - + — 
- 2 2n 



1 < m < n . 



(4.3) 



(4.4) 



and for every other value of q there are non-phase eigenvalues for some value of the rapidi- 
ties. 

This would lead to the result that the spectrum of the 1-folded model is always complex 
except for the isolated points 



arg(g) 



TT 



1 ^ 1 
- 1 ± - 

2 V n 



(4.5) 



4.2.2 Higher folded models based on 

We investigated the folded transfer matrices using numerical diagonalisation and found 
that the regions of real and complex spectrum are identical in every case to those of the 
1-folded models, independent of the folding number. 

It is obvious that the regions of complex spectrum contain those of the 1-folded model, 
since the spectrum of the folded model contains that of the 1-folded model, but it appears 
that there are no further constraints arising from the twisted sectors. 
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4.2.3 RSOS models based on 03 



The first observation to make is that since the spectrum of the RSOS models is a subset 
of the spectrum of the appropriate folded model, if the folded model has a real spectrum 
then so does the RSOS model. 

Next, since the RSOS spectrum is a subset of the folded spectrum, it is possible that the 
RSOS spectrum is real while that of the folded is complex. This is well known to be the case 
for the perturbations of the unitary minimal models A^p,p±i + 0i2 (note that 0i5 is never 
a relevant perturbation of a unitary minimal model), and we believe that this property is 
shared by many other models. 

As evidence for our assertion that the spectra of the RSOS models are subsets of those of 
the folded models, in figures H and |^ we plot the eigenvalues of the reduced two-particle 
transfer matrices for Kq-Kq scattering of the 10-folded model and of the associated RSOS 
models (We performed similar checks for the Kq-Ki and Kq-Kq-Kq transfer matrices with 
equally convincing results). 

These figures also show that for this particular choice of rapidity difference, there are 
regions of arg(g)/7r for which the folded transfer matrix has non-phase eigenvalues which 
are included in the regions for which we believe that the folded transfer matrices have 
non-phase eigenvalues for some values of rapidity. We also see that some of the RSOS 
restrictions do manage to omit these non-phase eigenvalues, while others do not. We 
report more fully on our findings in the later sections. 

4.2.4 The 0i2 perturbations with 71/7 > 1 

These theories have at least one higher kink Ki and we find that the most stringent 
restrictions already arise from the two-particle transfer matrix involving Kq and Ki. 

The theories with p' = ±1 mod p are of the 'unitary type': their S'-matrices (and transfer 
matrices) are proportional to those of the unitary models A4p^p+i + 0i2 and M.p^p+i + 02i, 
and since the scattering itself is manifestly unitary the eigenvalues are all phases. Indeed 
these S'-matrices do satisfy Hermitian analyticity which together with i?-matrix unitarity 
implies the S'-matrix unitarity equation for the two-particle S'-matrices. 

Our numerical calculations show that the transfer matrices of the theories A^py + 0i2 
where p' = {p ± l)/2 mod p also have phase eigenvalues only, although we do not know 
any explanation for this fact yet. 

Every other theory has non-phase eigenvalues in the Kq-Ki transfer matrix; conversely, 
the models described above still have phase eigenvalues when we consider the three-particle 
transfer matrices for all combinations of Kq and Ki. 

We have not tested these models beyond the three-particle transfer matrix. In addition, it 
is possible that particles arising through the S'-matrix bootstrap will also introduce non- 
phase eigenvalues in the transfer matrices; since the bootstrap has not been completed for 
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0.2 0.4 0.6 0.8 



1 



Figure 2: Comparison of the eigenvalues of the matrix part of the two-particle transfer matrices 

(2) 

for the RSOS model (A^io,m + (t>i2) and the 10-folded 03 ATFT The horizontal axis is arg(g)/7r 
restricted to the fundamental domain between and 1, and the vertical axis is the argument of 
the eigenvalues. The phase eigenvalues of the folded models are shown in green/thin solid lines, 
the non-phase eigenvalues in pink/thick solid lines, and the RSOS models (which are only defined 
for discrete values of q and happen in this case all to be pure phases) are shown as blobs. The 
relative rapidity of the particles was chosen 6 = 5. 

the majority of these models, we cannot say anything more about them. 
4.2.5 The perturbations with 1/2 < ir/'y < 1 

These theories are more usually known as the '02i' perturbations. They have no higher 
kinks in their spectra and at most one breather, and the S'-matrix bootstrap closes on these 
particles. The only possible violation of reality of the spectrum could be introduced by the 
fundamental kink S'-matrix, since the kink-breather and breather-breather S-matrices are 
pure phases. 

We show a selection of our results (up to folding number 50 and particle number 5) in 
figure |. 

The small (red) points are values which we find definitely to be 'non-real'. The remaining 
large (black) points are values which appear to have real spectra up to and including 
five-particle states. The solid (red, blue and purple) lines are the first three pairs of real 
series discussed in point 3 below and the dashed (green) lines are models for which a TBA 
equation for the ground state is known (see below). 
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Figure 3: Comparison of the eigenvalues of the matrix part of the two-particle transfer matrices 
for the RSOS model (A^5,m + </'i5) and the 10-folded af'' ATFT. The graph is labelled as for figure 
except that there are now also non-phase eigenvalues of the RSOS model which are shown as 
black blobs. The relative rapidity of the particles was chosen 6 = 5. 

We find a complicated pattern of results, with increasingly more models being ruled out 
by higher and higher particle number transfer matrices. Since we have only investigated 
transfer matrices involving at most 5 particles, our results are somewhat sketchy but they 
show the following patterns: 

1. All theories for which the folded models have real spectra are obviously real after 
RSOS restriction, that is those with |arg(g)/7r| = (1 ± l/n)/2. For the range of 7 
allowed, this gives only the models Ai4m,2m+i + 4>i2 = A^2m+i,4m + 02i and the models 

M + (f)l2=M 

m+l,2m+l + 021- 

2. The perturbations of the unitary minimal models have manifestly unitary S-matrices 
and real spectra, corresponding to the models Aip+i^p + (pu = Aip^p+i + 02i- 

3. It appears that several other infinite series of models may also have real spectra, at 
least this appears to be the case from our numerical tests of the transfer matrices 
up to 5 particles. These models form sequences which tend to the 'real' theories 
described in 1. above. In fact the first example of such a sequence are the unitary 
models which tend towards the value 71/7 = 1, which is the first 'real' model with 
n = 1. 

The next examples are the series ^^4^-1-1 3^-1-1 + tending towards arg(g)/7r = 3/4, 
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Figure 4: Numerical results for the eigenvalues of the transfer matrices of Mrs + 4'12- The 
heavy (black) points are models for which the multi-particle spectrum appears to be real up to 5 
particles, and the light (red) points those which are non-real, with [vr/7] = [s/r] plotted vertically 
against r= (folding number). See text for further details 



M3k±i,2k±i + (t>i2 tending towards arg(g)/7r = 2/3, Al8fc±3,5fe±2 + (t>i2 tending towards 
arg(g)/7r = 5/8, and so on... 

These three pairs of series and the unitary models are shown on figure ^ as solid lines. 

Some of these models have been considered before in various contexts. The first of these 
theories, 7^3,5 + (/>2i, was considered by G. Mussardo in [^], where he noted that its 



spectrum was real despite the fact that the S'-matrix was 'non-unitary'. 

Since then, (at least) three series of models have been conjectured to have a ground state 



described by real TBA equation^, these being M.m+i,2m+i + 02i IH, H^, -^2m+i,4m + 02i 
28|| , and M.2m+i,4:m-2 + 021 [ill- These are massive counterparts of the the |/|=1, 2 and 
4 massless flows considered in |^ with n/'j = 1/2 + |J|/(2r) and shown in figure || as 
dashed (green) lines. However, as is well known, simply knowing the TBA equations for 
the ground state does not determine the excited state spectrum uniquely, (a particular 
example of this effect being the 'type IT ambiguity noted in 0) and the scaling function 
can describe the ground state of a genuine unitary model and also a 'non-real' non-unitary 
model. So, while TBA equations are a useful alternative way to study the spectrum of 
a theory, it is not easy to deduce reality properties from the TBA equation for a ground 
state. 



'We thank R. Tateo for pointing this out to us. 



19 



4.2.6 The 0i5 perturbations with 7r/7 > 1 



These theories again have at least one higher kink Ki and again we find that the most 
stringent restrictions already arise from the two-particle transfer matrix involving Ko and 

Our numerical analysis shows that only two series of models have the possibility to be real. 
These are the theories with p' = 2p±l mod 4p for which the S'-matrix is Hermitian ana- 
lytic which guarantees the reality of the spectrum at the level of the Bethe-Yang equations, 
and the theories with p' = ±1 mod 4p for which the reality of the spectrum is mysterious. 

4.2.7 The 0i5 perturbations with n/'j < 1 

Again, these theories have no higher kinks in their spectra, and at most one breather, and 
the S'-matrix bootstrap closes on these particles. The only possible violation of reality 
of the spectrum could be introduced by the fundamental kink S-matrix, since the kink- 
breather and breather-breather S-matrices are pure phases. 
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Figure 5: Numerical results for the eigenvalues of the transfer matrices of Airs + </'i5- The 
heavy (black) points are models for which the multi-particle spectrum appears to be real up to 4 
particles, and the light (red) points those which are non-real, with 7r/7=s/4r plotted vertically 
against r=(folding number)/2. See text for further details 
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We find an even more complicated pattern of results than for the 0i2 perturbations, once 
again with increasingly more models being ruled out by higher and higher particle number 
transfer matrices. In this case the size of the transfer matrices is much larger, and we have 
only investigated transfer matrices involving at most 4 particles. Our results are even more 
sketchy but they show the following patterns: 

1. All theories for which the folded models have real spectra are still obviously real after 
RSOS restriction, that is those with |arg(g)/7r| = (1 ± l/n)/2. For the range of 7 
allowed, this gives only the models A^p,2p+i + 0i5 and M.2p+i,ip+i + 0i5- 

2. Again, it appears that several other infinite series of models may also have real 
spectra, at least this appears to be the case from our numerical tests of the transfer 
matrices up to 4 particles. These models form sequences which tend to the 'real' 
theories described in 1. above. 

The first example are the theories, J^k,ik-i + 0i5 tending towards arg(g)/7r = 1 and 
the series M.k,zk±i + 0i5 tending towards arg(g)/7r = 3/4. 

At the next level, we find that up to four particles, there are four infinite series 
tending to arg(g)/7r = 2/3, namely -Msfcii.sfcis + 0i5 and A^6fc±i,i6fc±2 + 0i5- Since 
the existence of this second series is in some sense a new phenomenon, we checked 
that it survives our numerical tests at the five particle transfer matrix level, but we 
have no opinion whether or not it will survive at all higher particle numbers. 

Again, three infinite series of these models have been considered before in various con- 
texts |2^. These series are exactly those which are related to the infinite series of 02i 
perturbations by interchanging the two terms in the ZMS potential as described in e.g. 0: 

M + 012 ^ M 2m+l,4m.+4 + 015 ; 

M 

4m,2m+l 

+ 012 ^ Mrn,2m+1 + 015 , (4.6) 

-M4m-2,2m+l +012 ^ M 2m— l,4m+2 + 015 • 

These series are shown as dashed (green) lines in figure ||. The final pair of theories is 
particularly interesting as being an example of a 'type 11' pair in that these are different 
theories but share exactly the same ground state scaling function and have a common 
sub-sector of multi-particle states with identical finite size energy corrections. 

5 Conclusions 

We have investigated the finite volume spectra of theories based on a^^^and Og^^ imaginary 
coupled affine Toda field theories. Our main results can be summarised as follows: 
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1. All theories (original, folded, RSOS) based on al have real spectra. To show this, we 
used the fact that the finite volume spectrum of the RSOS theories can be obtained 
as suitable projections of folded theories, which are in turn manifestly unitary QFTs 
related to sine-Gordon theory. We also presented new evidence for this (previously 
known) correspondence between the folded and RSOS models based on the transfer 
matrix method. 

(2) 

2. We conjectured that a similar correspondence exists between folded models and 
RSOS models based on , and supported this by numerical diagonalisation of their 
transfer matrices. 

3. We presented substantial evidence that unrestricted (both the folded and the original) 

(2) 

02 models have complex finite volume spectra in general, perhaps with the exception 
of a few special values of the coupling constant. 

(2) 

4. Similarly, it seems that RSOS theories based on 03 have complex spectra in gen- 
eral, with the exception of some special sequences asymptotically approaching the 
special values of the coupling for which the unrestricted models may have real spec- 
trum. These sequences in particular include the perturbations of the unitary minimal 
models, for which we know that the spectrum is real. 



There are quite a few open questions remaining. First of all, the NLIE description must be 

extended to excited states of 02 related models; this could provide us conclusive evidence 

(2) 

on whether or not the spectra of the unrestricted 03 theories are real for the special values 
of the coupling constants for which the transfer matrix calculations found no violation of 
reality. On the other hand, our results that show that the spectrum is complex in general 
must be reproduced by such an extension of the NLIE. 

Second, the projection of the folded spectrum to the RSOS spectrum must be explicitly 
realized. Together with the extension of the NLIE to excited states, this would open the 
way to (a) a systematic description of the spectra of the RSOS models; (b) determining 
whether the sequences for which we found no violation of reality with our methods really 
have real spectra. 

It seems appropriate to mention that another interesting problem is the relation of the 03 
model and '(/)2i' perturbations to the corresponding g-state Potts models. See |^ for TBA 
equations for the ground states of these models and for a discussion of the particle spectra 
in the g-state Potts models. 

(2) 

Third, as we mentioned in the text, the full spectrum of related theories is not yet 
closed in full generality. An example is Ads^^ + 0i2, for which the closed spectrum is not 
known up to date. One could attempt to close the S matrix bootstrap (we know of some 

attempts that failed, and we ourselves tried unsuccessfully); but equally well, it could 

(2) 

possibly be found by extending the NLIE to describe the excited state spectrum of 03 
related theories. 
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Appendix 



A Scattering amplitudes for Al^y + ^1^2 

The scattering amplitudes for the fundamental kinks in A^py + $1,2 were written down 
by Smirnov in ||12[. Here we briefly recall his result and give explicit expressions for the 
necessary q-6j symbols as it seems there are none available in the literature that are free 
of misprints and errors. 







222 Jvciax 2 Jvnax 2 

Figure 6: Adjacency rules for kinks in $12 perturbations for j^ax £ Z 

The allowed vacuum sequences {ji, . . . , j„} have the following adjacency rules: 
jfc+i = 1 if jfc = 

jk+i e {max(0, jfc - 1) • • -min (jfc + 1, p - 3 - jfc)} if jfc 7^ 



(A.l) 



and are of two types: either all jk are integers or all G Z + 1/2. The maximum value for 
jk is jmax = {p — 2)/2. The adjacency rules for these two sequences are shown on Figure 
H for the case when jmax is integer; for the case jmax G Z + 1/2 two very similar graphs 
result. Let us introduce the notation 

-'^ ^2tt^\ , . 2 



q = exp 



P 



X 



exp 



where ^ 



Tcp 
3 2p' — p 



(A.2) 
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and the g'-numbers 



qn_q n 

q-q-K 



Then the 2-kink scattering amplitudes take the form 

cd .^ _ {<r- l){<l' + l) , ,^fl_,\ J-UcM-c.-c.) fled 



where 



So{^) = ±^ (^sinh I - TTi) sinh | (^i? - ^ x 



oo 



sin ki) sinh ^ cosh (f - |) 

exp -2i / ^ \^^J—dk 

k cosh ^ sinh ^ 







Ca = a(a + 1) 

and the q-6j symbols take the form 

1 a-2 a-l1 _ f 1 a + 2 a + 1 
1 a a — 1 \ \ 1 a a + 1 



= 1 

la+la + ll [1 a a \ _ V [2a + 4] [2a] 



1 a a j ^ [la+la + lj^ [2a + 2] 
1 a a + 1 1 [2] 



1 a a + 1 [2a + l][2a + 2] 

1 a a - 1 1 [2] 



1 a a-1 [2a][2a + l] 

la+1 a 1 [1 a a + 1 1 _ 7[2a + 4][2c 
1 a a+lj^~\l a + 1 a j ^ ~ [2a + 2] 

la a — 11 |laa+l| _ ^[2a- l][2a + 3] 



laa + lj^ l^laa — IJ^ [2a + 1] 



la a 1 _ f 1 a a + 1 1 _ [2] / [2a + 3] 
1 a a+1 1 a a J ^ ~ ~ [2a + 2] y [2a + 1] 

1 a a1_ria + lal_ [2] 
1 a + 1 a J ^ \ 1 a a j ^ [2a + 2] 
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la a \ ( 1 a a-1 \ [2] / [2a - 1] 



1 a a - 1 j ^ [la a j ^ [2a] V [2a + 1] 

1 a al Jla — lal [2] 

1 a — 1 a j ^ \ 1 a aj^ [2a] 

1 a a 1 [2a- l][2a + 3] - 1 



1 a a J [2a] [2a + 2] 



(A.7) 



We remark that the square roots in (|A.7|) carry a sign ambiguity which must be resolved 



by an appropriate choice of the branches of the square root function. One must consider 
the square roots of each g-numbers separately and fix a sign for the expressions 



[n], n = l,2,... (A.8) 

and then use the selected representative consistently in all formulas. Such a choice of 
branch is necessary in order for the amplitudes to satisfy the Yang-Baxter equation, R- 
matrix unitarity and appropriate crossing relations. 



B Scattering amplitudes for Al^y + ^1^5 

In this appendix we present a corrected version of the S-matrix of $1,5 perturbations writ- 
ten down in the paper by one of the authors. The original formulas have misprints 
and some of them have the wrong normalisation factors. The correct ones can be obtained 
by imposing crossing symmetry and RU on the amplitudes. The choice of the normalisa- 
tion factors amount to defining the scalar product of the multi-kink states correctly (i.e. 
satisfying the constraints imposed by the quantum group symmetry). 

In $15 perturbations, the allowed vacuum sequences are composed of highest weights of 
the group Ug4,{sl{2)) where 

^ = (^) • (^-'^ 
They are labelled by jk = 0, 1/2, ... , jmax with jmax = (p — 2)/2. The adjacency rules for a 
sequence {ji, . . . , j„} are \jk+i ~ jk\ = or 1/2 (Figure |^). The zero difference corresponds 
to a neutral kink, while the nonzero to charged ones. Let us introduce the notation 

y = exp [ -t} ] , 4 



The scattering amplitudes for charged kinks are (with a sign corrected with respect to ||^) 



^ ^ q g+V V[2« + l]4[2c+l]4, 
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2"^ jrciax 2 jmax 

Figure 7: Adjacency rules for kinks in $15 perturbations 



where 5*0(^9) is the function defined in ( |A.5|) . For the sign ambiguities occurring as a result 



of square roots of g-numbers see the remark at the end of Appendix |^. 

The remaining amplitudes include neutral kinks. The neutral kink-neutral kink scattering, 
neutral kink-charged kink forward scattering and neutral kink-charged kink reflection were 
correct in the original paper and are the following 



eaa _ Q^y"^ + y'^(f -Q^- (fy"^ + y'^- o^^y"^ + yW - y'^q^ e i q^ 

^aa - 

^ba _ {y' + <fW-^) o^.. 

^ab — 51 Oo[tf , 

St = J-t^2M±Ils„W. (BA) 

The amplitudes describing two charged kinks turning into two neutral ones were incorrectly 
normalised (they include a nontrivial Clebsh-Gordan in the unrestricted theory Q). The 
correct form is 

C.aa _ ■ Caa {q^ - l){y^ - I) . 

^^'~'^a '^y ^^-^^ 

For the reverse process the amplitude is 

C.ab -Cab 

Saa = ^ 2 

Caa q^y 

where the normalisation factors Cab are 



Cab=i "i("^)'vlSk . (B.7) 

a2 a = b 

«! 2 are constants which are left free by the constraints of RU and crossing. The Cab appear 
in the crossing relations, which take the form 

S:t{^) = '^S'^lizTT - ^) . (B.8) 

^bc 
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C Transfer matrices and the Bethe-Yang equations 



C.l RSOS transfer matrices 



When we put the theory on a cylindrical space-time with spatial volume L, the allowed 
sequences of vacua (see appendices and P for 0i2 and 0i5 perturbations, respectively, 
and ||Tl| for the 0i3 case) are further restricted by the condition ai = a^+i, where 
is the number of particles. We take all particles to be point-like and ignore all vacuum 
polarisation contributions. Let us define the following (generalised) transfer matrix 



^2 

^1 ^2 



N 



which translates to 



'N, 



N 



a\a2...aisi 



(C.l) 



with the identification cat+i = Oi , hN+i = hi. In the following we shall not always write 
down the matrix indices explicitly. From the Yang-Baxter equation, it is straightforward 
to prove that these transfer matrices form a commuting family 



[T {§\§i, ...,^n),T (^9Vi, • • • , ^n)] = (C.2) 
We define the following specialised transfer matrices 

Apart from a phase coming from a plane wave factor of rapidity "i^fc in the wave function, 
this gives the monodromy corresponding to taking the fcth kink around the spatial circle 
multiplied by a factor (—1)"^. The total phase must equal (—1)^ depending on the statistics 
of the particle k {F = 1 for fermions F = for bosons). Thus we obtain the so-called 
Bethe-Yang equations |[T8| : 

exp (zm,i?sinh^,) ^ (^i, ^2, • • • , ^NttlX = (-l)^+^vl/^ife-''iv (c.4) 

where \[/"i"2 - <iiv jg ^j^g wave function amplitude, defined by the decomposition of the state 
|\E') as follows 

|V^)= J2 ^'''"''\Ka,aA^l)---Ka^aA^N)) . (C.5) 
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The energy and the momentum of the state (relative to the vacuum) are given by 

N N 

E = rrik cosh , P = sinh {^^ ■ (C-6) 

k=l k=l 

Due to the commutation relation ( |C.2| ), the equations ( |C.4D for the vector \1/ are compat- 
ible and can be reduced to scalar equations by simultaneously diagonalising the transfer 
matrices. Let us denote the eigenvalues of T . . . , by A*^*) • • • , ^n) with the 

corresponding eigenvectors ip^^^ i'&i, ■ ■ ■ ,''~^n) {s is just an index enumerating the eigenval- 
ues and the eigenvectors can be chosen independent of -d due to the commutativity (|C.2|) ). 
Then the solutions of the Bethe-Yang equations ( |U.4| ) are given by 

^aia2...a^ = (^1, • • • , ^ j^y^"^-"^ (C.7) 

where the rapidities solve the scalar Bethe Ansatz equations 

exp {tmkRsmh^k) A« {^k\A, ...,^n) = (-1)^+' (C.8) 



C.2 Folded transfer matrices 

For folded models the vacua are labelled by an integer a modulo the folding number k. 
The allowed sequences satisfy 

aj_|_i = tti + Q mod k (C9) 

where Q are the possible topological charges of the solitons 

/ ±1, sine - Gordon , . 

^ ~ \ +1, 0, -1 ZMS ^^-^"^ 

In finite volume we require Oi = Qn+i mod k for periodic boundary conditions. Then the 
transfer matrix has the same form as in equation ( |C.1| ) except that now the matrix S is 
constructed from the scattering matrix S of the sine-Gordon/ZMS model in the following 
way: 

where Qab denotes the charge of the soliton connecting the vacua a and b (see Figure H). 
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Q 



ad 



dc 




ab 



Q 



be 



Figure 8: Vacuum labels and topological charges for two-particle ^-matrices in folded 
models 



References 



[1] C. Mudry, C. Chamon and X.-G. Wen, Nud.Phys. B 466 (1996) 383-443, 
|cond-mat/9509054 . 



N. Read and H. Saleur, Nucl. Phys. B 613 (2001) 409-444, |hep-th/0106T24 



M. J. Bhaseen, J.-S. Caux, I. I. Kogan and A. M. Tsvelik, Nucl. Phys. B 618 (2001) 



465-499, |cond-mat/0012240 



[2 
[3 
[4 
[5 
[6 

[7 
[8 
[9 



G. Takacs and G. Watts, Nucl. Phys. B547 (1999) 538-568, |hep-th/9810006 . 
P. Dorey and R. Tateo, Nucl. Phys. B 571 (2000) 583-606, |hep-th/9910T0^ . 
G. Takacs, Nucl. Phys. B 489 (1997) 532-556, |hep-th/9604098 



J.L. Cardy and G. Mussardo, Phys. Lett. B225 (1989) 275-278. 

H.G.Kausch, G.Takacs and G.M.T.Watts, Nucl. Phys. B 489 (1997) 557-579, 
|hep-th/ 9605104 . 



D.I. Olive, Nuovo Cimento 26 (1962) 73. 

J.L. Miramontes, Phys. Lett. B455 (1999) 231-238, |hep-th/990rTi5 



T.R. Klassen and E. Melzer, Nucl. Phys. B 338 (1990) 485-528. 



[10] Z. Bajnok, L. Palla, G. Takacs and F. Wagner, Nucl. Phys. B 587 (2000) 585-618, 
|hep-th/0004181 . 



[11] N. Reshetikhin and F.A. Smirnov, Commun. Math. Phys. 131 (1990) 157-178. 

[12] F.A. Smirnov, Int. J. Mod. Phys. A6 (1991) 1407-1428. 

[13] A. Kliimper and F.A. Pearce, J. Stat. Phys. 64 (1991) 13; 

A. Kliimper, M. Batchelor and P.A. Pearce, J. Phys. A24 (1991) 3111. 

[14] C. Destri and H.J. de Vega, Phys. Rev. Lett. 69 (1992) 2313; Nucl. Phys. B 438 
(1995) 413-454, [hep-th/9407lT7[ . 



29 



[15] D. Fioravanti, A. Mariottini, E. Quattrini and F. Ravanini, Phys. Lett. B 390 
(1997) 243-251, |hep-th/960809l| . 



[16 
[17 

[18 
[19 

[20 

[21 

[22 
[23 
[24 

[25 
[26 

[27 

[28 

[29 

[30 



C. Destri and H.J. de Vega, Nud. Phys. B 504 (1997) 621-664, [hep-th/970lT07 . 

G. Feverati, F. Ravanini and G. Takacs, Phys. Lett. B430 (1998) 264-273, 
[hep-th/9803T04[ ; Nucl. Phys. B 540 (1999) 543-586, [hep-th/9805TT7| ; Phys. Lett. 



B444 (1998) 442-450, |hep-th/9807T60i 



T.R. Klassen and E. Melzer, Nucl. Phys. B 382 (1992) 441-485, |hep-th/9202034 



RE. Dorey and R. Tateo, Nucl. Phys. B 482 (1996) 639-659, [hep-th/9706T40| ; Nucl. 
Phys. B 515 (1998) 575-623, |hep-th/9607167[ . 

Al. B. Zamolodchikov, Nucl. Phys. B 432 (1994) 427-456; Phys. Lett. B335 (1994) 
436-443. 

G. Feverati, F. Ravanini and G. Takacs, Nucl. Phys. B 570 (2000) 615-643, 
[hep-th/990903l . 



R Zinn- Justin, J. Phys. A31 (1998) 6747-6770, [hep-th/9712222 



S. Coleman and H.J. Thun Commun. Math. Phys. 61 (1978) 31. 

V.B. Lidskii, U.S.S.R. Comput. Math, and Math. Phys. 6 (1966) 73-85, 
J. Moro, J.V. Burke and M.L. Overton, SIAM Journal on Matrix Analysis and 
Applications 18 (1997) 793-817. 

G. Mussardo, Int. J. Mod. Phys. A7 (1992) 5027-5044. 

F. Ravanini, M. Stanishkov and R. Tateo Int. J. Mod. Phys All (1996) 677-698, 
|hep-th/9411085| . 



M.J. Martins, Phys. Rev. Lett.69 (1992) 2461-2464, [tiep-tli/9205024| ; 
Nucl. Phys. 5394 (1993) 339-355. 

E. Melzer, Supersymmetric Analogs of the Gordon- Andrews Identities, and Related 
TEA Systems, [hep-th/9412l54[ . 

RE. Dorey, C. Dunning and R. Tateo, Nucl. Phys. B 578 (2000) 699-727, 
[hep-th/OOOllSS . 



P.E. Dorey, A. Pocklington and R. Tateo, Integrable aspects of the scaling g-state 
Potts models I: bound states and bootstrap closure and Integrable aspects of the 
scaling g-state Potts models II: finite size effects (to appear). 



30 



